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基于 移 不 变 字 典 学 习 和 稀 琉 编码 的 滚动 轴承 故障 识别 算法 
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摘 要 : 针对 现 有 旋转 机 械 故 障 识别 算法 过 度 依 赖 专 家 先 验 知识 的 问题 ， 提 出 了 一 种 基于 移 不 变 字典 学 习 和 稀疏 编码 


(SIDL-SC) 的 自 适 应 故障 识别 算法 。 首 先 ， on 下 的 振动 信号 进行 分 段 和 平滑 预 处 理 以 降低 数据 处 理 复 
， 接 着 将 加 入 自 适 应 惩罚 因子 的 移 不 变 字 典 学 习 算 法 用 于 提取 不 同 故障 状态 下 的 移 不 变 基 函数 ; 然后 ， 利 用 高 效 


Se ge ni 的 重 构 ; 最 后 ， 以 重 构 残 差 作为 
对 该 信号 故障 状态 识别 的 判断 依据 。 滚 动 轴承 振动 数据 库 和 实测 航空 发 动机 振动 信号 的 实验 结果 表明 ， 该 算法 相 较 于 
现 有 算法 具有 更 高 的 故障 识别 准确 率 ， 在 实际 中 具有 较 强 的 可 行 性 。 
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Fault recognition algorithm for rolling bearings based on 
shift invariant dictionary learning and sparse coding 
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Abstract: According to current algorithms for rotating machines largely depending on expert prior knowledge, the paper 
proposed an adaptive fault recognition algorithm based on shift invariant dictionary learning and sparse coding. Firstly, it 
segmented and smoothed vibration signals to decrease the complexity. Then, it used shift invariant dictionary learning with 
adaptive penalty factor to learn shift invariant bases in different fault states. After that, it used an efficient sparse coefficient 
solver called Feature Sign Search for reconstructing signal to be recognized. Lastly, residual was an evidence to determining 
fault state the signal belonging to. In the experiments of rolling bearing datasets and vibration signals of real aero-engine 
demonstrate its higher accuracy than up-to-date algorithms and feasibility for practical applications. 
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承 故 障 下 的 特征 频率 和 包 络 谱 ， 通 过 构建 线性 分 类 器 实现 不 同 
故障 的 识别 。 但 是 基于 频谱 的 分 析 方 法 依然 依赖 于 人 工 经 验 知 


引言 


对 机 械 设 备 进行 状态 监控 和 故障 诊断 对 于 提高 生产 效率 、 识 ， 缺 乏 对 数据 自身 的 有 效 挖掘 。 王 录 雁 等 人 四 提出 了 
预防 事故 发 生 具 有 重要 的 实际 意义 。 在 机 械 故 障 诊断 中 ， 对 监 。 EMD (experience mode decomposition ) 和 灰色 关联 理论 的 轴承 
控 设 备 有 效 提取 故障 特征 是 实现 故障 模式 识别 的 重要 一 环 。 然 。 故障 诊断 方法 , 利用 EMD 将 振动 信号 进行 内 豪 模 态 分 量 (IMF ) 
而 ， 机 械 设 备 多 样 的 工 况 和 复杂 的 机 械 结构 给 特征 提取 带 来 了 分解, 建立 各 IMF 能 量 分 布 与 轴承 故障 状态 的 灰色 关联 模型 实 


难度 ， 基 于 时 域 和 频 域 的 特征 提取 算法 依赖 于 人 工 特征 选择 ， 岗 故 障 诊断 。 但 是 EMD 自身 存在 模 态 混 闭 现象 ， 对 于 信号 不 


基于 


~ 


FE 需 要 相关 领域 内 的 专家 知识 才能 较 好 的 完成 特征 提取 工作 。” 同 状态 下 的 能 量 分 布 可 能 会 产生 干扰 。 蔡 剑 华 国 等 提出 了 


主 介 


高 阶 统计 量 的 滚动 轴承 故障 诊断 方法 ， 该 方法 的 实质 是 将 信号 
息 


一 直 以 来 ， 研 究 人 员 针 对 机 械 振动 故障 特征 提取 做 了 大 量 ”的 高 阶 谱 恢复 为 功率 谱 ， 利 用 功率 谱 提取 信和 号 的 故障 特征 
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信 ， 


的 工作 。 其 中 ，Harmouche 等 人 DI 利 用 频谱 分 析 法 研究 不 同 各 实现 故障 诊断 。 但 是 高 阶 统计 量 存在 受 奇异 值 影 响 较 大 的 缺陷 
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导致 算法 的 结果 并 不 稳定 。 
近年 来 基于 字典 学 习 的 特征 提取 算法 成 为 机 器 学 习 领 域 的 


研究 热点 ， 其 核心 思想 是 通过 字典 学 习 在 原始 数据 中 构建 稀 下 C 
变换 域 , 使 待 识别 信号 在 该 变换 域 下 的 投影 具有 稀疏 性 ,Wright 


等 人 中 提出 了 基于 稀 玻 表示 分 类 ( sparse representation 
classification, SRC) 的 人 脸 识别 算法 , 首次 在 稀疏 表示 理论 与 
模式 识别 之 间 建 立 了 纽带 ， 使 基于 稀疏 表示 的 特征 能 够 作为 模 
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曲 建 岭 ， 等 : 基于 移 不 变 字 典 学 习 和 稀 朴 编码 的 滚动 轴承 故障 识 另 
标 函 数 
Ce 
其 中 : 稀 琉 惩罚 项 8 用 于 控制 稀 玻 化 程度 。 通 常 1 设 为 固定 值 ， 


但 由 于 训练 样本 长 度 的 不 同 ， 固 定 值 往往 并 不 能 达到 最 优 解 ， 
因此 本 文 考虑 对 其 进行 改进 以 适应 在 不 同 训练 样本 长 度 下 p 
能 具有 较 好 的 控制 稀 疏 化 效果 。 令 8 =09.L， 其 中 工 为 训练 样 


式 识别 的 特征 向 量 实现 对 样本 模式 的 准确 识别 ， 这 给 故障 模式 
识别 带 来 了 新 的 思路 。 

字典 的 构建 方式 通常 分 为 两 种 : 基于 数学 分 析 的 字典 构建 
方法 和 自 适 应 字典 构建 方法 。 基 于 数学 分 析 的 字典 构建 方法 针 
对 原始 数据 的 先 验 知识 实现 稀疏 化 ， 常 以 正 交 变换 基 的 形式 实 
现 过 完备 字典 的 构建 。 振动 数据 常 以 离散 余弦 变换 中 (DFT) 和 
傅 里 叶 变 换 四 (FFT) 作为 正 交 变换 基 。 使 用 上 述 这 些 正 交 变换 
基 虽 然 具 有 相对 完善 的 算法 框架 ， 但 缺乏 灵活 性 ， 不 能 保证 对 
应 的 数据 类 型 在 该 变换 基 下 是 最 稀 疏 的 表示 ， 某 些 信 号 特别 是 
时 频 变 化 范围 较 广 的 信号 在 上 述 变换 基 下 仍 不 稀疏 或 不 够 稀 玻 ， 
从 而 影响 对 数据 的 重 构 精 度 四 。 基 于 自 适 应 的 字典 构建 方法 则 
以 原始 数据 为 训练 样本 ， 通 过 K- 奇 异 值 分 解 00 (K-SVD ) 等 算 
法 实现 过 完备 字典 的 构造 ， 往 往 采用 这 种 方法 构造 的 字典 对 于 
原始 数据 的 表达 更 具有 稀疏 性 。 从 模式 识别 角度 来 看 ， 越 稀 玻 
的 特征 往往 对 应 着 更 佳 的 分 类 效果 。 
对 于 旋转 机 械 设 备 而 言 ,同一 特征 的 出 现 具有 循环 往复 性 ， 
因此 本 文 考虑 将 移 不 变 字 典 学 习 算 法 引入 旋转 机 械 故障 识别 中 ， 
以 期 能 从 原始 数据 中 捕捉 到 具有 时 移 不 变性 质 的 本 质 特征 构成 
基 函 数 ， 将 这 些 特征 平移 来 构建 过 完备 字典 ， 实 现 对 机 械 设备 
故障 状态 的 有 效 识别 。 


1 ”SIDL-SC 故障 识别 算法 理论 


1.1 移 不 变 特征 自学 习 算 法 

2006 年 ，Smith0 首 先 在 《Nature》 杂 志 上 提出 了 移 不 变 稀 
跑 编码 学 习 算 法 (shift invariant sparse coding，SISC) 并 成 功 应 用 
于 声 信号 特征 提取 。 在 形式 上 ,SISC 理论 认为 同一 特征 模式 在 
间 序 列 上 的 移动 是 一 个 独立 事件 ， 学 习 到 的 基 函 数 具 有 相对 
时 间 的 移 不 变性 。 在 构成 学 习 字 典 的 原子 上 , SISC 学 习 到 的 基 
函数 是 信号 中 反复 出 现 的 本 质 特征 ， 原子 体现 了 信号 
a le 因此 ,SISC 良好 的 移 不 变 特性 使 其 已 在 图 
像 处 理 03、 信 和 号 处 理 03] 等 领域 得 到 了 应 用 。 

SIsc 使 用 个 包含 M 个 基 函 数 的 字典 [gd ,gq,，,. 
对 于 任意 输入 信号 y e ;* 均 可 由 这 些 林 函数 和 相应 
积 和 来 表示 ， 基 本 数学 模型 为 


至 


每 个 字 


对 ,| si QxM 
系数 的 卷 


M 
= Dd,*X,,+e (1) 
m=] 


一 


其 中 :“* ”表示 卷 积 算 子 ， 用 于 表示 基 函 数 在 所 有 可 能 位 置 
的 移动 。 系 数 x ，e ;具有 稀疏 性 ， se 通常 表示 高 斯 白 噪声 。 
各 参数 满足 以 下 关系 : O<NW，P<N 且 P+O-1=N。 定 义 目 


本 长 度 ，g 定义 为 稀 玻 化 尺度 〈 通 常 0<g<0.1)， 由 此 实现 对 
惩 如 项 的 改进 。 

字典 学 习 阶 段 若 同时 优化 基 函 数 o 和 系数 x， 则 该 优化 
问题 为 非 凸 的 ， 无 法 得 到 稳定 解 。 因 此 ， 


间 


借鉴 稀 玻 编码 中 交替 
多 化 的 思想 , 依次 更 新 o 和 X ，， 反复 迭代 直到 目标 函数 收敛 。 
包 化 过 程 分 为 系数 求解 和 字典 学 习 两 部 分 049。 


在 系数 求解 部 分 ， 固 定 字典 姜 求 解 稀疏 系数 x ， 则 式 (1) 模 
型 被 分 解 为 ;=1,N 共 N 个 独立 的 优化 问题 , 针对 每 个 输入 信 
号 y 求解 得 到 的 系数 都 独立 于 其 他 的 输入 ， 因 此 考虑 单个 输入 
y 的 系数 求解 优化 问题 。 


M 2 
minly, -20, 2 


本 文采 用 高 效 的 Batch-OMP 算法 E15], 
佳 的 匹配 原子 x ， 其 表达 式 为 


{x'} 全 argmin ly, -dx 
[ * 


在 字典 万 中 寻找 最 


tol < vx (4) 


在 字典 学 习 部 分 ， 固 定 系数 x 求解 字 
问题 转换 为 乙 约束 优化 问题 : 


和 矩阵 Dp， 则 原 优化 


win -> > *X,, | sz lo <c,1<j<M (5) 
直接 求解 上 式 计 算 量 巨大 ， 本 文 在 字典 学 习 阶段 借鉴 时 频 


域 转换 的 思想 ， 将 时 域 计 算 量 较 大 的 卷 积 运算 转换 到 频 域 ， 
频 域 的 乘积 运算 代替 时 域 的 卷 积 运算 ， 使 运算 复杂 度 由 O(N?) 
降低 为 O(Nlog N) ， 提 高 了 运算 效率 。 由 Parseval 定理 可 知 ， 
式 (5) 可 转换 为 如 下 优化 问题 19: 


> 


N 


min2 必 


c 为 常数 ，g 为 基 函 


sto se (0) 


其 中 数 dl 的 离散 传 里 叶 变换 ， 相 应 地 ， 


和 和 了 分 别 为 x 和 y 的 离散 传 里 叶 变 换 。 使 用 拉 格 朗 日 乘 子 
法 求解 该 问题 ， 则 该 问题 被 分 解 为 二 次 项 求 和 的 形式 ， 对 于 每 
一 个 频 域 上 的 频率 / ， 构 造 拉 格 朗 日 函数 : 


za=2 
其 中 ，4ei "为 对 侦 变量 ， 单 位 向 量 /。,; "， 并 且 


EF wt th 
MM Xa aa) 
Of) A 
,A=diag(N) ,X=| CD ED 
de! 
yw 


YY 一 XN 


+(O0Y AG®]- cr 0O) 


ey 
a 
~ 
dl 二 » yy = 
(1) 
dy 
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六 


L(9,7) 是 复数 函数 ， 通 过 复数 公式 | 王 |ReCz)| +|Im(z)| 


将 复数 函数 转换 为 实 函数 求解 ， 将 Rec@) 和 fm 计算 目标 函 


数 0m =argminL(9, 力 ,优化 参数 4 可 利用 牛顿 法 求解 ,代入 式 


G@) 即 实现 对 基 的 迭代 更 新 。 
day @ 


1.2 基于 特征 符号 搜索 的 系数 求解 算法 

考虑 到 Batch-OMP 算法 需要 预先 设 定 稀疏 度 等 参数 , 这 在 
一 定 程度 上 增加 了 算法 的 计算 量 ， 如 果 参 数 选择 不 当 很 可 能 收 
敛 速度 较 慢 甚至 难以 收 勾 。 考 虑 到 待 识 别 样本 的 稀 玻 度 一 般 未 
知 ， 因 此 本 文 在 测试 样本 识别 阶段 引入 特征 符号 搜索 算法 
(Feature Sign Search，FS) 求 解 待 识别 信号 在 不 同 基 函 数 上 的 重 
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因此 考虑 将 重 构 残 差 作 为 故障 识别 的 判断 准则 。 利 用 FS 算法 
求解 待 识 别 信号 在 各 个 状态 移 不 变 基 函数 下 的 稀 朴 系数 
由 式 (11) 求 解 基 函数 及 其 对 应 的 稀疏 系数 的 重 构 


bb 


六 多 和 


少 = ya sx (11) 


寻 此 ，SIDL-SC 算法 的 实现 分 为 数据 预 处 理 、 移 不 变 字 
学 习 、 重 构 系 数 求解 和 计算 残 差 四 部 分 ， 有 具体 实现 步骤 如 图 1 
所 示 。 


测试 样本 类 别 


测试 样本 


fs 系数 求解 


1 > 重 析 入 


移 不 变 字 FS 系数 求人 重 构 信号 2 


典 2 


故障 状态 2 一 | 分 段 | “| 移 不 


构 系 数 。 FS 算法 由 Lee 提出 09,， 其 核心 思想 是 通过 对 系数 符号 
的 猜测 将 系数 分 量 分 为 可 行 集 (@ ) 和 0 和 集 (@ ),， 在 0 集中 
不 断 搜索 梯度 变化 最 大 的 分 量 并 将 其 对 应 的 下 标 i 放 入 可 行 集 
中 直到 算法 收 全 ,算法 具体 步骤 如 下 : 

a) 初 始 化 x=0，0=0， 其 中 6 ={-1,0,1} 表示 系数 分 量 x 
的 符号 。 

b) 在 系数 为 零 的 分 
放 入 可 行 集中 。 


中 搜索 梯度 变化 最 大 的 分 量 及 其 下 标 


赂 


ugme Lv 00 -ou (9) 


oy 


若 >y， 则 9= 1; 着 


oy -4x| <-y， 则 
OX 


oly-AY4 
OX, 


0 =1。 
0) 选 取 和 矩阵 A 的 子 矩 阵 A， 及 中 只 含有 与 可 行 集中 i 相对 
应 的 列 向 量 , 同样 ,选取 与 ;相应 的 子 系数 向 量 和 ,计算 无 


约束 优化 方程 min (9 =| -AR| +7y6 福 ， 求 得 最 优 解 


多 ,=(AA (ANy -yO" /2) (10) 


在 % 和 ,之 间 进行 线 搜索 ， 找 到 所 有 符号 改变 的 系数 向 
量 , 并 将 它们 的 目标 函数 与 在 &_ 的 目标 函数 进行 比较 , 将 又、 
修正 为 使 得 目标 函数 值 最 小 的 系数 向 量 。 将 中 系数 分 量 为 夫 


的 ; 从 可 行 集 移 除 ， 并 更 新 符号 


合 g。 
9 着 wx #0， 都 有 时 4 ,ysignoc) -0 成立, 风 转 
X 


条 件 (b)， 否 则 转 6); 


四 若 Yx =0， 都 有 时 二 人 |<，， 则 得 到 最 优 解 x， 否 
X 


则 返回 b)。 

1.3 ”基于 重 构 残 差 的 故障 识别 算法 

理论 上 来 说 ， 与 待 识别 信号 的 真实 状态 越 接 近 ， 则 由 基 函 
数 及 其 求解 的 稀 玻 系数 重 构 的 信号 与 待 识别 信号 的 误差 越 小 ， 


人 65 系数 求 多 重 构 信号 N 


Y 


图 1 本 文 故 障 识别 算法 流程 图 


故障 状态 N 一 > 


2 实验 


2.1 实验 设置 

本 文选 取 美 国 凯 斯 西 储 大 学 的 开放 轴承 振动 数据 库 进 行 实 
仿 11， 使 用 电 火 花 加 工 技术 分 别 在 轴承 外 圈 、 内 圈 和 滚动 体 上 
设置 单 点 故障 , 故障 直径 分 别 为 0.18 mm, 0.36 mm, 0.54 mm。 
本 文 将 不 同 故 障 数 据 分 为 两 组 ， 如 表 1 所 示 。 其 中 实验 1 使 用 
数据 集 A 中 的 样本 进行 训练 和 测试 , 实验 2 中 采用 实验 1 的 字 
息 学 习 到 的 基 函 数 对 1~2 hp 负载 下 的 样本 进行 测试 , 以 检测 本 
文 算法 的 泛 化 能 力 。 


表 1 实验 数据 集 描述 


数据 集 ”负载 hp ”故障 类 型 ”故障 尺寸 /mm ”类 别 标记 
A/B 0/1~3 正常 0 1 
圈 故 障 0.18 2 
圈 故 障 0.36 3 
团 故 障 0.54 4 
外 圈 故 障 0.18 5 
外 圈 故 障 0.36 6 
外 圈 故 障 0.54 7 
滚动 体 故障 0.18 8 
滚动 体 故 障 0.36 9 
滚动 体 故障 0.54 10 


2.2 SIDL-SC 模型 的 建立 

参数 的 选择 影响 算法 的 计算 效率 和 故障 识别 的 准确 度 ， 医 
此 有 必要 根据 训练 样本 的 实际 情况 筛选 出 各 个 参数 的 最 优 值 作 
为 实现 较 高 状态 识别 率 的 保证 , 本 节 考 察 对 于 建立 SIDL-SC 模 
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型 起 重要 作用 的 几 个 参数 。 识别 率 。 对 于 滚动 轴承 来 说 ， 学 习 到 的 基 函 数 应 该 包含 轴承 中 
2.2.1 稀疏 化 尺度 至 少 一 个 完整 的 冲击 ， 从 实际 测量 来 看 ， 滚 动 轴承 的 两 个 相 邻 

在 移 不 变 字 典 学 习 中 ， 本 文采 用 自 适 应 稀疏 尺度 参数 g 探 波峰 的 间隔 平均 在 76 左右。 同样 以 上 一 节 中 的 实验 样本 为 例 ， 


制 稀 疏 系 数 在 目标 函数 中 的 比重 ,防止 系数 过 大 而 基 函 数 过 小 。 ”不 同 基 函 数 长 度 下 的 样本 识别 率 如 图 2(b) 所 示 ， 从 图 中 可 以 看 
实际 表现 为 较 大 的 0 值 具有 较 快 的 收敛 速度 但 是 较 高 的 重 构 误 出， 基 函 数 在 70~90 能 够 达到 较 高 的 样本 识别 率 。 而 基 函 数 长 
差 。 图 2(a) 以 实验 集 A 中 第 2 类 内 圈 18mm 故障 样本 为 例 ， 比 。 度 较 短 时 , 由 于 每 个 移 不 变 基 函数 不 能 学 习 到 完整 的 信号 冲击 ， 
较 了 不 同 g 值 下 移 不 变 字典 学 习 重 构 的 原始 信号 的 残 差 。 从 图 ”导致 正确 率 不 能 达到 最 佳 。 而 基 函 数 长 度 较 长 可 能 会 带 来 信息 
中 可 以 看 出 当 稀 琉 太 度 0 较 低 时 ， 有 具有 较 低 的 重 构 误差 ， 随 着 。 元 余 , 不 利于 后 期 使 用 FS 进行 重 构 。 因 此 ,本 文选 择 基 函数 长 
稀疏 尺度 的 增 大 ， 重 构 误 差 呈现 显著 增 大 的 趋势 。 但 是 从 计算 。 上 度 为 80。 

效率 来 说 ， 0 较 小 带 来 了 每 次 迭代 需要 的 计算 量 较 大 ， 从 而 导 ”2.2.3 训练 样本 长 度 的 选择 

致 收敛 速度 较 慢 。 因 此 从 重 构 误差 和 计算 时 间 效 率 两 方面 来 均 在 多 次 实验 中 发 现 不 同 的 训练 样本 长 度 对 于 故障 识别 的 作 
衡 ， 选 择 稀疏 尺度 为 0.025 是 合适 的 。 很 大 。 以 0 hp 负载 下 第 1、2、5、8 类 样本 为 例 进行 训练 和 
测试 。 其 中 每 类 故障 状态 随机 选择 100 个 样本 用 于 测试 ， 实 验 


~ 


1 | 结果 如 图 3 所 示 。 从 图 中 可 以 明显 看 出 , 当 训练 样本 长 度 为 128 

让 可 和 256 时 ， 在 第 8 类 故障 的 重 构 中 ， 第 2、5、8 类 学 习 到 的 基 

/A | 函数 都 能 够 以 较 低 的 残 差 对 其 重 构 ， 具 体 表现 在 这 三 类 基 函 数 
Eo 的 重 构 误 差 出 现 了 混 基 ， 这 给 样本 的 正确 识别 带 来 了 困难 。 而 
| ce | 输入 样本 为 512 和 1024 则 较 好 的 避免 了 混 炙 的 现象 ， 通 过 对 
人 是 比 原始 信号 发 现 ， 训 练 样本 长 度 为 128 和 256 分 别 包含 1 和 3 


0 50 80 7 80 9 10 10 120 
基 函 数 长 并 


个 完整 冲击 , 而 长 度 为 512 和 1024 包含 有 更 多 的 完整 冲击 , 从 
基 函 数 自学 习 的 角度 来 看 ， 较 多 的 冲击 能 使 基 函 数学 习 到 的 特 
征 更 加 显著 ， 更 有 利于 后 期 的 重 构 分 类 。 但 训练 样本 长 度 的 增 
2.2.2 基 函 数 长 度 的 选择 加 会 带 来 每 批 次 学 习 时 间 的 增长 ， 因 此 选择 训练 样本 长 度 为 
在 移 不 变 字 典 学 习 中 ， 基 函数 的 长 度 的 选择 影响 对 样本 的 512。 


(a) 稀 玻 尺度 对 重 构 误 差 的 影响 (b) 基 函数 长 度 与 识别 准确 率 的 关系 
图 2 稀 琉 尺度 和 基 函 数 长 度 对 重 构 误 差 的 影响 


ER 


InputLength = 128 


InputLength = 256 


Reconstruction Eror 


Reconstruction Error 


Reconstruction Error 
四 > 四 [3 己 总 


图 3 训练 样本 长 度 与 重 构 误 差 的 关系 

2.3 ”实验 结果 表 2 SIDL-SC 各 参数 值 

经 过 上 一 节 的 实验 讨论 , 最 后 确定 了 SIDL-SC 算法 中 各 主 参数 数值 
要 参数 值 如 表 2 所 示 。 图 4 展示 了 部 分 学 习 到 的 移 不 变 基 函数 ， ee 
在 实验 一 中 分 别 在 数据 集 A 中 每 类 状态 随机 选择 200 个 样本 用 ei 
于 测试 ， 实 验 结果 如 表 3 所 示 。 Po 

从 表 3 可 以 看 出 ， 本文 所 提出 的 方法 在 数据 集 A 中 的 表现 a 
良好 , 大 部 分 样本 都 可 以 达到 100% 的 正确 识别 率 , 只 有 部 分 样 样本 重 受 率 50% 
本 存在 1~2 个 识别 错误 ， 随 着 数据 集 样本 量 的 增 大 ， 本 文 算法 de 
可 以 取得 更 高 的 正确 识别 率 。 


201805.00216v1 


chinaXiv 


2 人 


Chi 
等 : 基于 移 不 变 字典 学 习 和 的 滚动 轴承 故障 识 


录用 稿 曲 建 岭 ， 
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SAA 
| 
100 200 300 400 500 600 700 800 900 1000 0 10 2 3 40 的 
. = HENESE 
o 8) 
图 4 各 状态 典型 振动 信号 及 其 对 应 的 移 不 变 基 函数 
(a(b). 正 常 状 态 ;，c(qd). 内 圈 故 障 ; (9. 外 圈 故 障 ;，g(D). 滚 动 体 故 障 ) 
人 和 稀 玻 编码 框架 内 实现 对 轴承 故障 的 识别 ， 其 未 探 讨 参数 选择 
样本 类 别 ”识别 正确 率 (%) ”样本 类 别 “识别 正确 率 (%) 问题 ， ra ke 7 人 该 元 余 
， 和 站 号 的 标准 差 、 偏 度 、 哨 度 等 时 、 括 吉村 征 送 入 郑 积 神经 网 络 
， 和 ， 证 分 类 器 ,利用 LeNet5 模型 实现 对 故障 的 识别 。 将 以 上 两 种 方法 
的 和 本 文人 法 者 在 0hp 下 的 样本 进行 训练 分 别 在 数据 集中 各 
从 实验 结果 可 以 看 出 ， 采 用 传统 时 频 域 特征 用 于 分 类 的 CHEN 
2.4 与 传统 方法 对 比 实 验 的 算法 具有 相对 较 低 的 识别 率 , 而 LIU 的 算法 虽然 在 个 别 样本 


为 了 体现 SIDL-SC 算法 在 轴承 故障 识别 中 的 优越 性 , 本 文 
与 其 他 部 分 学 者 算法 进行 了 对 比 实验 。LIU08 同 样 在 字典 学 习 
1 et | 
3 : 
HEE 样本 类 别 0 Rt 8 1 
(al hp 下 识别 率 (b)2 hp 下 识别 率 


2.5 ”实测 航空 发 动机 实验 


图 5 


不 同 负载 下 各 算法 样本 识别 率 


上 的 识别 率 较 高 ， 但 总 体 不 如 本 文 算法 的 稳定 性 高 。 


(c)3 hp 识别 率 


为 验证 本 文 算法 在 实际 装备 中 的 可 行 性 ， 将 该 算法 应 用 于 #9 es 
实测 某 型 航空 发 动机 故障 识别 。 实 验 数据 采集 于 该 型 直升机 涡 
轴 发 动机 多 台 次 、 多 次 台 架 试车 时 所 测 得 的 四 种 状态 下 发 动机 ~ A 
振动 信号 ， 即 正常 状态 、 转 子 不 平衡 故障 、 转 子 松动 故障 和 碰 0 0 
摩 故障 ， 采 用 等 弧度 采样 方式 ， 每 周期 采集 128 个 数据 点 ， 典 VA AA 
型 振动 信号 如 图 6 所 示 。 0 

利用 SIDL-SC 算法 对 不 同 状态 下 的 振动 信号 进行 移 不 变 = < 

参数 设置 如 下 训练 样本 长 度 为 256， 基 函数 长 度 2 
为 50， 稀疏 化 尺度 为 0.01， ee 其 余 参数 设置 与 人 
表 3 相同 。 学 习 到 的 移 不 变 基 函 数 如 图 7 所 示 。 
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图 7 不 同 状 态 下 学 习 到 的 移 不 变 基 函 数 〈a.j 


从 图 中 可 以 看 出 正常 状态 下 的 移 不 变 基 函数 波形 比较 平缓， 
故障 状态 下 学 习 到 的 移 不 变 基 函数 均 不 同 程度 的 含有 其 各 自 独 
特 的 特征 。 对 不 同 状态 进行 识别 ， 每 组 实验 重复 进行 100 次 ， 
取 识 别 准 确 率 作为 评价 指标 ， 实 验 结果 如 表 4 所 示 ， 并 与 文献 
[18] 和 文献 [19] 中 的 方法 进行 对 比 。 实 验 结果 表明 ， 由 于 发 动机 
转子 不 平衡 故障 和 转子 松动 故障 与 正常 状态 的 波形 上 存在 一 些 
相似 之 处 , 导致 算法 的 识别 率 低 于 CWRU 数据 集 , 但 相 较 于 上 
述 两 种 算法 ， 本 文 算法 仍然 能 保持 较 高 的 故障 识别 率 ， 进 一 步 
验证 了 本 文 算法 的 可 行 性 。 


表 4 航空 发 动机 实验 识别 结果 


识别 率 (%) ”正常 状态 ”转子 不 平衡 ”转子 松动 ” 碰 摩 
文献 [18] 方 法 80 78 74 86 
文献 [19] 方 法 76 73 74 80 
本 文 算法 85 82 81 93 


3 ”结束 语 


SIDL-SC 是 一 种 结合 移 不 变 字典 学 习 和 稀 玻 编 码 的 振动 信 


号 故障 识别 算法 ， 尤 其 是 针对 多 工 况 条 件 下 的 故障 识别 有 着 较 
高 的 正确 识别 率 。 该 方法 以 原始 振动 数据 作为 训练 样本 ， 通 过 
特征 自学 习 捕 获 反 映 不 同 状 态 信 号 本 质 特 征 的 移 不 变 基 函数 ， 


以 最 小 重 构 残 差 为 准则 判别 故障 的 类 别 。 同 时 对 移 不 变 字典 委 
习 中 的 惩罚 因子 进行 了 优化 ， 以 稀 玻 尺度 作为 前 者 的 因 变 量 ， 
有 效 解 决 了 不 同 训练 样本 长 度 下 惩罚 因子 的 设 定 问 题 。 以 
CWRU 振动 数据 库 作为 研究 对 象 , 本 文 详细 探讨 了 不 同 参数 选 
择 对 实验 结果 的 影响 ， 送 选 出 最 优 参数 作为 最 终 的 实验 参考 变 
量 。 同 时 ， 将 本 文 算法 应 用 于 实测 航空 发 动机 振动 信号 ， 并 与 
其 他 部 分 学 者 的 算法 进行 实验 对 比 ， 结 果 表 明 ， 本 文 算法 具有 
更 强 的 泛 化 能 力 ， 可 以 更 好 的 应 用 于 实际 振动 数据 的 处 理 。 
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